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Gaps in the spectrum of the Neumann Lapla- 
cian generated by a system of periodically dis- 
tributed trap 

Andrii Khrabustovskyi 1 2 3 , Evgeni Khruslov 2 

■ Abstract. The article deals with a convergence of the spectrum of the Neumann Laplacian in a periodic 

unbounded domain Q. £ depending on a small parameter e > 0. The domain has the form D. E = R" \ S £ , 
where S £ is an eZ' ! -periodic family of trap-like screens. We prove that for an arbitrarily large L the spectrum 
has just one gap in [0, L] when e small enough, moreover when e — > this gap converges to some interval 
whose edges can be controlled by a suitable choice of geometry of the screens. An application to the theory 
of 2D-photonic crystals is discussed. 

Keywords: periodic domain, Neumann Laplacian, spectrum, gaps 

^ Introduction 

' 

It is well-known (see, e.g., 0]-[3]|) that the spectrum of self-adjoint periodic differential operators 
has a band structure, i.e. it is a union of compact intervals called bands. The neighbouring bands 
may overlap, otherwise we have a gap in the spectrum (i.e. an open interval that does not belong to 
\ the spectrum but its ends belong to it). In general the existence of spectral gaps is not guaranteed. 

For applications it is interesting to construct the operators with non-void spectral gaps since 
their presence is important for the description of wave processes which are governed by differen- 
tial operators under consideration. Namely, if the wave frequency belongs to a gap, then the cor- 
responding wave cannot propagate in the medium without attenuation. This feature is a dominant 
requirement for so-called photonic crystals which are materials with periodic dielectric structure 
attracting much attention in recent years (see, e.g., EH6]). 

In the present work we derive the effect of opening of spectral gaps for the Laplace operator in 
R" in > 2) perforated by a family of periodically distributed traps on which we pose the Neumann 
boundary conditions. The traps are made from infinitely thin screens (see Fig. [T|). in the case n = 2 
this operator describes the propagation of the //-polarized electro-magnetic waves in the dielectric 
medium containing a system of perfectly conducting trap-like screens (see the remark in Section 
©. 

We describe the problem and main result more precisely. Let e > be a small parameter. Let 
S E = IJS^bea union of periodically distributed screens S E in R" (n > 2). Each screen S ? is an 

(n - l)-dimensional surface obtained by removing of a small spherical hole from the boundary of 
a n-dimensional cube. It is supposed that the distance between the screens is equal to s, the length 
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of their edges is equal to be, while the radius of the holes is equal to de 1 ^- if n > 2 and e l/dE ~ if 
n = 2. Here d e (0, oo), 6 (0, 1) are constants. 



Fig. 1. The system of screens S E 

By J\ E we denote the Neumann Laplacian in the domain W \ S E . Our goal is to describe the 
behaviour of its spectrum cr(Jl E ) as s — » 0. 

The main result of this work is as follows (see Theorem ll.ll ). For an arbitrarily large L the 
spectrum cr(jR E ) has just one gap in [0, L] when s is small enough. When e — > this gap converges 
to some interval (cr,fj.) depending in simple manner on the coefficients d and b. Moreover (see 
Corollary 11.11 ) with a suitable choice of d and b this interval can be made equal to an arbitrary 
preassigned interval in (0, oo). 

The possibility of opening of spectral gaps by means of a periodic perforation was also inves- 
tigated in [7]. Here the authors studied the spectrum of the Neumann Laplacian in K 2 perforated 
by Z 2 -periodic family of circular holes. It was proved that the gaps open up when the diameter 
d of holes is close enough to the distance between their centres (which is equal to 1). However, 
the structure of the spectrum in []7] differs essentially from that one in the present paper. Namely, 
when d — > 1 the spectrum converges (uniformly on compact intervals) to a sequence of points. 

Various examples of scalar periodic elliptic operators in the entire space with periodic coeffi- 
cients were presented in [|8l-fT8l. In these works spectral gaps are the result of high contrast in 
(some of) the coefficients of the operator. 

The outline of the paper is as follows. In Section \T\ we describe precisely the operator Jl E and 
formulate the main result of the paper (Theorem II .It describing the behaviour of cr(J{ £ ) as e — » 0. 
Theorem [Tj] is proved in Section |2] Finally, in Section |3]on a formal level of rigour we discuss the 
applications to the theory of 2D photonic crystals. 

1. Setting of the problem and the main result 

Let n e N \ { 1 } and let s > 0. We introduce the following sets: 

• B = {x = (xi, . . . , x„) 6 R" : -b/2 < x\ < b/2, Vi'}, where b e (0, 1) is a constant. 

• D E = [x 6 OB : |jc — *°| < d E ), where x° = (0, 0, ... , 0, b/2) and d E is defined by the fol- 
lowing formula: 



de nl , 



s exp 



ds 2 



n>2, 
n = 2. 



(1) 



Here d > is a constant. It is supposed that s is small enough so that d E < b/2 



• S £ = dB\D £ 
For i 6 Z" we set: 



S £ = s(S £ + i) (2) 

and (see Fig. [T]) 



Q e = R" \ 



VieZ" > 



(3) 



Now we define precisely the Neumann Laplacian in Of. We denote by tj e [u, v] the sesquilinear 
form in L 2 (Q e ) which is defined by the formula 

rf[u, v] = J (Vw, Vv) dx (4) 

, _ v-i <9v 

and the definitional domain dom(?7 £ ) = H (Q E ). Here (Vw, Vv) = > . The form rj £ [u, v] is 

j-* dx k dx k 



densely defined closed and positive. Then (see, e.g., |fT91 Chapter 6, Theorem 2.1]) there exists the 
unique self-adjoint and positive operator 3\ £ associated with the form rf\ i.e. 

(Ji £ u, v) Lim = 7j £ [u, v], Vm e domC?l £ ), Vv e dom(?7 e ) (5) 

It follows from © that Jl £ u = -Au in the generalized sense. Using a standard regularity theory 
(see, e.g, EOl Chapter 5]) it is easy to show that each u e dom(J?I) belongs to Hj oc (Q. £ ), furthermore 
du 

— = for any smooth Y c dQ. £ . 
on r 

We denote by cr(Jl £ ) the spectrum of Jl £ . To describe the behaviour of cr(Jl £ ) as s —* we need 
some additional notations. 

In the case n>2we denote by cap(T') the capacity of the disc 

T = {x = (x u ...,x n ) e R" : \x\ < 1, x n = 0} 

Recall (see, e.g, [21]) that it is defined by 



cap(7) = inf f \Vw\ 2 dx 
w J 

where the infimum is taken over smooth and compactly supported in K" functions equal to 1 on T. 
We set 

°- = \ , Abn n = ^-r (6) 

\w> n = 2 > 

It is clear that cr < /j. 

The behaviour of <t(3\ £ ) as s —* is described by the following theorem. 

Theorem 1.1. Let L be an arbitrary number satisfying L > [x. Then the spectrum of the operator 
J\ £ has the following structure in [0, L] when s is small enough: 

(rC?l £ )n[0,L] = [0,L]\(cr £ ,p £ ) (7) 
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where the interval (<x £ ,/i e ) satisfies 



limcr e = cr, lim// = u (8) 



Corollary 1.1. For an arbitrary interval (o~,p) c (0, oo) there is a family {Q. £ } £ of periodic un- 
bounded domains in W 1 such that for an arbitrary number L satisfying L > p the spectrum of the 
corresponding Neumann Laplacian 3\ E has just one gap in [0, L] when s is small enough and this 
gap converges to the interval (cr, p.) as s — > 0. 

Proof It is easy to see that the map (d, b) i-» (cr,p) is one-to-one and maps (0, oo) x (0, 1) onto 
{(<x,yu) 6 K 2 : < <x < //}. The inverse map is given by 

r^4 T (l-^-.)(cap(T))-., "> 2 - b= ^7P (9) 
|2cr7r (1 — o~p~ ), n = 2, v H 

Then the domain Q e considered in Theorem [Tj] with d and b being defined by formula © satisfies 
the requirements of the corollary. □ 

Remark 1.1. It will be easily seen from the proof of Theorem II . 1 I that the main result remains valid 
for an arbitrary open domain B which is compactly supported in the unit cube (-1/2, 1/2)" and 
whose boundary contains an open flat subset on which we choose the point x°. In this case the 
coefficients cr and p are defined again by formula © but with |Z?| instead of b" (here by | • | we 
denote the volume of the domain). 

Remark 1.2. One can guess that in order to open up m > 1 gaps we have to place m screens S £ , 
5 2,- • • ,S £ m in the cube (-1/2, 1/2)". However the proof of this conjecture is more complicated 
comparing with the case m = 1 . We prove it in our next work. Below we only announce the result 
for the case m = 2. 

Let B\ and B 2 be arbitrary open cuboids which are compactly supported in (-1/2, 1/2)". On 
dB i and dB 2 we choose the points x 1 and x 2 correspondingly. We suppose that these points do not 
belong to the edges of cuboids. Let D E and D e 2 be open balls with the radii d\ and d E and the centres 
at x 1 and x 2 correspondingly. Here dp j = 1,2 are defined by formula (Q3 but with dj instead of d 
(dj > are constants). For / e Z", j = 1, 2 we set Sf- = s(S £ + i) and finally 



Of = R" \ 



u »•« 

\ieZ"J=\,2 



By Jl E we denote the Neumann Laplacian in Q £ . 

We introduce the numbers cry, j = 1,2 by formula © with dj and \Bj\ instead of d and b" 
correspondingly. We suppose that dj are such that the inequality cr\ < cr 2 holds. Finally we define 
the numbers pj by the formula 

= J (Pi + Pi + CT\ + CT 2 + (-I) 7 V(Pl +P2 + O-1+ Cr 2 ) 2 - 4(p 1 CT 2 + PiCT 2 + O^O^j 

where pj = cr ; -|fi ; -|(l - \Bj\y 1 . It is not hard to check that o~\ < p\ < cr 2 < p 2 . 

Now, let L be an arbitrary number satisfying L > p 2 ■ Then cr(Jl E ) has just two gaps in [0, L] 
when e is small enough, moreover the edges of these gaps converge to the intervals (crj,pj) as 
s — > 0. Also it is easily to show that the points crj,pj can be controlled by a suitable choice of the 
numbers dj and the cuboids Bj. 
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2. Proof of Theorem I1.1I 

We present the proof of Theorem 11.11 for the case n > 3 only. For the case n = 2 the proof is 
repeated word-by- word with some small modifications. 

In what follows by C, C\... we denote generic constants that do not depend on e. 
By {u) B we denote the mean value of the function v(x) over the domain B, i.e. 

{u)b = j^j J u(x)dx 

B 

Recall that by \B\ we denote the volume of the domain B. 

2.1. Preliminaries. We introduce the following sets (see Fig. E): 

Y = {xeR": -1/2 < < 1/2, Vi}. 
Y £ = Y \5 e . 
F = Y\B. 

Let A £ be the Neumann Laplacian in e _1 Q £ . It is clear that 

a(Jl £ ) = £rV(A £ ) (10) 

A £ is an Z"-periodic operator, i.e. A £ commutes with the translations u(x) i-» w(jc + 0> 2 e F° r 
us it is more convenient to deal with the operator A £ since the external boundary of its period cell 
is fixed (it coincides with dY). 





D E 

/ 


















B 



Fig. 2. The period cell Y £ 

In view of the periodicity of A £ the analysis of the spectrum cr(A £ ) is reduced to the analysis 
of the spectrum of the Laplace operator on Y £ with the Neumann boundary conditions on S £ and 
so-called 6-periodic boundary conditions on dY. Namely, let 

T = {e = (9 u ...,e n )eC n : V£|0*| = 1} 

For 9 G T" we introduce the functional space H x e (Y £ ) consisting of functions from H l (Y £ ) that 
satisfy the following condition on dY: 

Vk=\,n: u(x + e^) = Qku(x) for x = (x\, x 2 , ■ ■ ■ , -1/2, . . . , x n ) (11) 

T 

fc-th place 

where e k = (0, 0, . . . , 1 , . . . , 0). 
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By rf- e we denote the sesquilenear form defined by formula © (with Y £ instead of Q.) and 
the definitional domain H l e (Y £ ). We define the operator A e ' £ as the operator acting in L 2 (Y £ ) and 
associated with the form rf ,£ , i.e. 

(A £ ' 9 u, v) L2(Y c } = r] £ ' 6 [u, v], Vw 6 dom(A £ ' 9 ), Vv e dom(rj £ ' d ) 

The functions from dom(A e ' £ ) satisfy the Neumann boundary conditions on S £ , condition (fTTI) on 
dY and the condition 

du du 

Vk=l,n: —(x + e k ) = 6 k —(x)forx = (xi,x 2 ,...,-l/2,...,x n ) (12) 

OX], OXk 

The operator A e ' £ has purely discrete spectrum. We denote by {/^' £ } the sequence of eigenval- 
ues of A e ' £ written in the increasing order and repeated according to their multiplicity. 

The Floquet-Bloch theory (see, e.g., Q3-T51) establishes the following relationship between the 
spectra of the operators A £ and A e ' £ : 

OO 

cr(A £ ) = \J L k , where L k = {J {A d k ' s } (13) 

k=l 6eT" 

The sets L k arc compact intervals. 

Also we need the Laplace operators on Y £ with the Neumann boundary conditions on S £ and 
either the Neumann or Dirichlet boundary conditions on dY = dY £ \S £ . Namely, we denote by 
rf' e {resp. rj D,£ ) the sesquilinear form in L 2 (Y £ ) defined by formula © (with Y £ instead of Q £ ) and 
the definitional domain H\Y £ ) {resp. H { (Y £ ) = [u £ H\Y £ ) : u = on dY £ \ S £ )). Then by A N > E 
(resp. A D ' £ ) we denote the operator associated with the form 77^ (resp. ?] D ' £ ), i.e. 

(A £ -*u, v) L2(n = t] £ '*[u, v], Vm e dom(A £ *), Vv e dom(7/ £ '*) 

where * is (resp. D). 

The spectra of the operators A N ' £ and A D ' £ axe purely discrete. We denote by (^f' e | (resp. 

{^k £ }k<Et? ^ e se£ l uence °f eigenvalues of A N,£ (resp. A D ' S ) written in the increasing order and 
repeated according to their multiplicity. 

From the min-max principle (see, e.g., [Q~l Chapter XIII]) and the enclosure H l (Y £ ) d H l e (Y £ ) d 
Hq(Y £ ) one can easily obtain the inequality 

Vfc e N, V# e T" : A N k £ < A e k £ < A D k £ (14) 

In this end of this subsection we introduce the operators which will be used in the description 
of the behaviour of A k , Af and A\ as e — » 0. By A^ (resp. A^, A^.) we denote the operator which 
acts in L 2 (F) and is defined by the operation A, the Neumann boundary conditions on dB and the 
Neumann (resp. Dirichlet, ^-periodic) boundary conditions on dY. By A B we denote the operator 
which acts in L 2 (B) and is defined by the operation A and the Neumann boundary conditions on 
dB. Finally, we introduce the operators A N , A D , A 6 which act in L 2 (F) © L 2 (B) and are defined by 
the following formulae: 

A N -- 



(1 






-(1 






ii 












As/' 






As) 



We denote by {^} km ( res P- { A k} km > l A k\ km ) the 

sequence of eigenvalues of A N (resp. A D , A 6 ) 
written in the increasing order and repeated according to their multiplicity. It is clear that 

Af = ^ = 0, A% > (15) 

A» = 0, ^ > (16) 

A? = 4 = 0, A| >0if = (1,1 1) (17) 

4 = 0, ^ > Oif 0*(1,1, ...,1) (18) 

2.2. Asymptotic behaviour of Dirichlet eigenvalues. We start from the description of the as- 
ymptotic behaviour of the eigenvalues of the operator A D - E as s — » 0. 



Lemma 2.1. For eac/? e N one /las 



Furthermore 



\xm^ = X D k (19) 

e— *u 



Af e ~ ere 2 as s -> (20) 

where cr is defined by formula (0). 

Proof. We start from the proof of (fT9l . It is based on the following abstract theorem. 

Theorem. (Iosifyan et al. 1122^) Let *H £ , "K be separable Hilbert spaces, let £ E : "H £ — > Ff £ , X° : 
— > < H° oe linear continuous operators, imX° c "V c Ff , where *V is a subspace in Ff . 
Suppose that the following conditions C\ - C4 hold: 

C\. The linear bounded operators R E : "K — > < H E exist such that ||^ £ /||L £ — > <c? | ] 1 1 an y 

£— >0 " 

/ e < V. Here g > is a constant. 

C 2 . Operators £ E ,£° are positive, compact and self-adjoint. The norms ||-£ £ HxcH £ ) cire bounded 
uniformly in s. 

C 3 . For any f e <V: \\£ E R E f - R E £°f\\w- -» 0. 

e— »0 

C4. For any sequence f E e *H £ such that sup ||/ £ ||<he < 00 the subsequence s' c e and w e *V 

£ 

exist such that \\£ E f E - R E w\\^. — > 0. 

£=£'—►0 

Then for any k e N 

£— >0 

where {jufjtl] {/^tl^Li are ^ e eigenvalues of the operators £ E and £ Q , which are renumbered in 
the increasing order and with account of their multiplicity. 

Let us apply this theorem. We set <H E = L 2 (Y E ), <H° = L 2 {F) © L 2 (B), £ E = (A D ' E + I)" 1 , 
£° = (A D + I)' 1 (here I is the identity operator), <V = Ff , the operator R E : ( H° — > , K £ is defined 
by the formula 

[R E n(y) = \{ F( fl' y£F D f = (fF,f B )e'H (21) 

Obviously conditions C\ (with q = 1) and C2 hold (namely, ||-£ E ||£(««) < 1). 
Let us verify condition C 3 . Let / = (f F , f B ) e -K . We set f E = R E f, v E = £ E f E . By v £ and v E B 
we denote the restrictions of v E onto F and 5 correspondingly. It is clear that 

\K\\ 2 HHF) + \k\\ 2 hHb) = \\v £ \\ 2 HI(Yn < nrwl^ = 211/11^ (22) 
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We denote 

= {we H\F) : w\ dY = o) 

Since Hq(F) © H l (B) is compactly embedded into 'H then due to estimate (1221) there is a subse- 
quence s' <z e and v F 6 H^(F), v B e H l (B) such that 

Vp — > v f weakly in H l (F) and strongly in L 2 (F) 

e=s'->0 (23) 
v £ — > v B weakly in Z/ 1 (5) and strongly in L 2 (B) 

One has the integral equality 

(Vv £ , Vw £ ) + v £ w £ - fw 



dy = 0, Vw £ e //^(y £ ) (24) 



ye 

We introduce the set 



W = |(w F , w fl ) 6 C°°(F) © C°°(5) : w F | f5y = 0, (supp(w F ) U supp(w B )) n {x } = 

(here as usual by supp(/) we denote the closure of the set {x : fix) £ 0}). Let w = iw F ,w B ) 
be an arbitrary function from W. We set w £ = R £ w. It follows from the definition of W that 
supp(w £ ) n Z) s = when e is small enough and therefore w E e C°°(Y £ ) n Hq(Y s ). 

We substitute w £ into (1241) and taking into account (1231 we pass to the limit in (1241) as s = s' —> 0. 
As a result we obtain 



(VVg, Vw B ) + V B W B - f B W B 



dy = (25) 



Since is a dense subspace of H^(F) © H l iB) then (1251) is valid for an arbitrary (w F ,w B ) e 
//q(F) © //'(-B)- It follows from ([25]) that -Afv F + v F = f F and -A B v B + v B = fs and consequently 

v = X /, where v = iy F , v B ) (26) 

We remark that v do not depend on the subsequence e' and therefore the whole sequence iy s F , v e B ) 
converges to (v F , v B ) as s — » 0. Condition C 3 follows directly from (|2TT >. (1231 and (|26T ). Obviously 
condition C4 was proved during the proof of C 3 . 

Thus the eigenvalues /i £ of the operator £ £ converges to the eigenvalues fi k of the operator X° as 
e -» 0. But e = OuP" 1 - 1, A* = Ou^)" 1 - 1 that implies ©. 

Now we focus on the proof of (1201 . Let vf' £ be the eigenfunction of A D £ that corresponds to the 
eigenvalue Af' E and satisfies 

llvf ' £ Hl 2( ^) = 1 (27) 

(vf ' E ) B > (28) 

One has the following inequalities 

||vf' £ ||i 2(F) < C\\Vv?>% (F) (29) 

\\v^-(vf' £ ) B \\l 2(B) <C\\Vvf'% 2{B) (30) 

\B\ ■ |<vf £ ) B | 2 < ||vf % (B) < 1 (31) 



Here the first one is the Friedrichs inequality, the second one is the Poincare inequality and the 
third one is the Cauchy inequality. Furthermore one has 

HVvf £ HL (F£) = (||vf' e ||^ + ||vf £ - (vf E ) B \\l 2(B) + \B\ ■ |<vf e > B | 2 ) (32) 

Below we will prove (see inequality (1561 )) that 

if £ < Cs 2 (33) 

Then it follows from <|37), (|29]>-([33]) that 

UVvf £ \\ 2 Ll(YE) < Cs 2 (34) 
||vf % (F) < Cs 2 (35) 
||vf £ - <vf £ > J 2 2(B) < Ce 2 (36) 
The last inequality can be specified. Namely, one has 

1 = llvf £ |lL (f) + llvf £ - (vf E ) B \\ 2 L2(B) + \B\ ■ |<vf S > B | 2 (37) 
Then in view of (|33T)-(|3"71) 

|i<vf e > s i 2 - i^r 1 ! < C^ 2 

Finally, taking into account ([281 ) we get: 



lim||vf' e -|5|- 1/2 ||L B1 = (38) 

ct a convenient ; 

following problem 



Now we construct a convenient approximation vf £ for the eigenf unction vf £ . We consider the 



A<A = in R" \ T (39) 
iff = 1 in 57/ (40) 
i/Kjc) = o(l) as |x| -> oo (41) 

Recall that I = jj £ R" : |jc| < 1, x„ = 0}, obviously <9r = T. It is well-known that this problem 
has the unique solution ifr(x) satisfying J \Vi//\ 2 dx < oo. Moreover it has the following properties: 

R»\r 

iff e C°°(R" \ T) (42) 

l/f(Xi ,X 2 ,..., Xn-u x n ) = if/(x\ , x 2 , . . . , X n -l, -x„) (43) 

cap(T) = J \Vi//\ 2 dx (44) 

R n \r 

The first two properties imply: 

— = in {x e W : x n = 0}\T (45) 

ox n 

dip dip 

dx n ,v„=+o dx n 

Furthermore the function tff(x) satisfies the estimate (see, e.g, Il23l Lemma 2.4]): 

|£>>(x)| < C\x\ 2 ~"- a for \x\ > 2, \a\ = 0, 1, 2 (47) 



= in T (46) 

x„=-0 
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We define the function vf' e by the formula 



vf' £ « 



1 



2^8] 
1 



•A 



x - x 



\x - X 



1 



1 



x - X 



\x - X 



x e F 
x£BUD E 



(48) 



where <p : 



i vrsr im 

4 R is a twice-continuously differentiable function such that 
ip(p) = 1 as p < 1/2 and <£>(p) = as p > 1, 
/ is an arbitrary constant satisfying 

< / < imin{l - b,b} 

Here we also suppose that s is small enough so that d E < 1/2. It is easy to see that the constructed 
function vf' £ (x) belongs to dom(A D - £ ) in view g0]), dH, g^), (O, (ill), (|50]>. 
Taking into account (1441) . (1471) we obtain: 



(49) 



(50) 



||Vvf £ || 2 , (FE) ~ 4- 1 cap(r)rf n - 2 |5|- 1 £ 2 = as 2 (s -> 0) 



,D,£||2 



<Ce 4 



(51) 
(52) 



Since vf' e = on dY and [vf £ - |flf 1/2 ] . nf = on dB \ [x : \x - x°\ < l) (here [. . . ] int means the 
value of the function when we approach dB from inside of B) we have the following Friedrichs 
inequalities 

llvf £ |] 2 2(F) < C||Vvf % 1(F) (53) 



\B\- XI2 \\l 2(B) < C||VvflH 2 



(54) 



It follows from (ED, (ED, d54j) that 

K- £ llL ( y,) ~ 1 (e^O) (55) 

Using the min-max principle (see, e.g., [Q] Chapter XIII]) and taking into account (ISTI) . (1551) we 
get 



,1 



IIVv 



D,s\\2 



< 



1 Hia(J*) 



cap(r)J"' 2 | J Sr 1 e 2 = cre z (e -> 0) 



1 || Z), £||2 



Now let us estimate the difference 



w £ = vf e -vf £ 



(56) 



(57) 



One has 



2{\\v?% (F) + \\yf'% (F) ) + 2(||vf- £ - IBI-^HU, + 111^51 1/2 - vf' % m ) 
and thus in view of ((35) . (1381) . (IBTI) . (1531) . (l54l) we get 



,£||2 



IL 2 (F 



< rv 2 



(58) 



Substituting the equality Vj e = v\ ' e + w e into (|56l) and integrating by parts we obtain 



||Vw e || 2 2(y£) <2(Avf e ,w £ ) L2(yE) + 



v iiv D ii 2 



IVv 



1 Ili2(7 s ) 
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and in view of (|5TT> . (l52l) . (l55l) . (l58l) we conclude that 



lim£f 2 ||Vw £ || 2 =0 (59) 

Finally using (l5TT) . (l59l) we obtain 

if £ ~ ||Vvf £ || 2 2(y£) ~ ^ (e ^ 0) (60) 

The lemma is proved. □ 

2.3. Asymptotic behaviour of Neumann eigenvalues. In this subsection we study the behaviour 
of the eigenvalues of the operator A N,e as s —* 0. 



Lemma 2.2. For each k e N one /zas 



limif = ^ (61) 



Furthermore 

A N 2 £ ~ us 1 ass^O (62) 

where fi is defined by formula ©. 

Proof. The proof of (16TT ) is similar to the proof of ([T91 , so we focus on the proof of (|62l . 
Let v^' £ be the eigenfunction of A N ' S that corresponds to A 2 ,£ and satisfies 

llvf Ww) = 1 (63) 

(vf ) B > (64) 
Since the eigenspace which corresponds to A N y ,£ = consists of constants then 

(vf , l)L 2( y s) = (65) 



and in view of the min-max principle we have 



HV< ( 



^ = 1^11^)= ™? -rr^ 1 (66) 



where //](7 £ ) = jv 6 ^(y e ) : (v £ , 1)^) = 0) 
Below we will prove (see inequality (1761 ) that 



A N 2 £ < Cs 2 



Then in the same way as we obtain the estimates for vf £ in Lemma [2TI we derive the estimates for 

vf: 

HVvf \\ 2 Li(¥S) < Cs 2 (67) 
llvf - (vf ) F ||L (F) + Hvf - (v N 2 £ ) B \\ 2 L2iB) < Cs 2 (68) 
Using d67J), d68]) we obtain from d63]» and d65J that 

|<vf > F | 2 |F| + |<vf > B | 2 |5| ~ 1 (s -> 0) 
<vf > F |F| + <vf>*|fl| = 
and therefore taking into account (l64l we get 



<vf> F ~ - Vl^l/IF, <vf ) B ~ j\F\/\B\ (s -> 0) (69) 



12 



We construct an approximation vf for the eigenfunction vf by the following formula 



N,£ 



vf (*) = 




1 



+ 



\F\ 2^/\B\\F\ \ d 



x — x- 



\\F\ 



x- X K 



x - X 



\x - X 



1 



xeF 



xeBUD £ 



(70) 



\B\ 2yf\B\\F\ 

Recall that iff is a solution of (T39~l)-(l4TI). <p : R — » R is a twice-continuously differentiable function 
satisfying ( |49l ), / is a constant satisfying (1501) . One can easily show (taking into account the equality 



Vvf ||i 2(ye) ~ 4- l \F\- l \B\- l c2Lp(T)d n - z s z = /ue z (s -» 0) 



-1 1 rri — 1 1 r>i — J 



n-2 2 



|Av 



\B\ + \F\ = 1) that vf e dom(A yv - £ ). 

In the same way as we obtain the estimates for vf ' £ in Lemma |2~T1 we get the estimates for vf : 

(71) 
(72) 

(73) 

(74) 
(75) 



2 "L 2 {Y- 



< Ce 4 



vf + 



^J\B\/\F\ 



L 2 (F) 



+ 



- (e -» 0) 



L 2 (B) 



< Ce 2 



||vf|| 2 ~l ( e ->0) 



Since vf - << ' e > e Hl(Y £ ) then it follows from d7D, dTJ), (ESJ that 



,lf = ||Vv 



< 



IIVv 



N,sn 2 
2 <>L 2 (Y £ ) 



2{YS) \\y N - 



/v N < E \ II 2 
(v 2 )rr" 



lis 2 (s -> 0) 



Now let us estimate the difference 



One has 



w £ = vf -vf 



(76) 



(77) 



vf + viwi 



L 2 (F) 



+ 



\\w% (YE) <2\ 



and thus in view of d68J), (|69l) . (1731) we get 



JV,s 



2 

L 2 (F) 



+ 2 



- j\F\/\B\ 



+ 

2 

L 2 (B*) 



\im\\w £ \\ 2 Li(Yn = 



£^0 



j\F\/\B\-v 



2 



Substituting the equality vf = vf + w £ into (1761) and integrating by parts we get 



||Vw £ || 2 2(y£) <2(Avf,w e ) L2(yE) + 



Me ...£\ 



IIVv 



N,su2 
2 M L 2 (F £ ) 



vf - (vf >, 



- HVvf || 



2 llz,2(y £ ) 



Vll'2 



and then in view of CD}, CZ2]), (f74l . (1751) . (1781) we conclude that 



lime- z ||Vw £ || 

s->0 



,e||2 _ q 



Finally using dTn) . (1791) we obtain 



A 



N,B 



HVvf llL ( y E) ~^ 2 (e->0) 



The lemma is proved. 



L 2 (B £ ) 



(78) 



(79) 

(80) 
□ 
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2.4. Asymptotic behaviour of ^-periodic eigenvalues. To complete the proof of Theorem 11.11 
we have study the behaviour of the spectrum of the operator A e ' E . 

Lemma 2.3. For each 6 e T" and ieN one has 

lim4 e = 4 (81) 

Furthermore 

A 2 E ~fis 2 i/0 = (l,l,...,l) (82) 
A e { £ ~ ae 2 i/0* (1,1,..., 1) (83) 

Proof. The proof of (18TI) is similar to the proof of (fT9l ). 

The proof of (1821) is similar to the proof of (l62l . Namely, we approximate the eigenfunction 



N,s 



v e 2 e of A ft£ that corresponds to A e 7 ,£ and satisfies Hv^Hz^y*) = 1, {v 6 2 ' £ )b > by the function v 2 
(170b (since v^' £ is constant in the vicinity of dY then it satisfies (fTT) - ([T2|) with = (1, 1)). The 
justification of the asymptotic equality 

4 £ ~ HVvf ||2 2(ye) , = (1,...,1) 

is completely similar to the proof of (ISUT l. 

Now, we focus on the proof of (18"31 . Let Vj' £ be the eigenfunction of A e ' e that corresponds to A 6 ^ 
and satisfies ||Vj' £ || L2 (y E ) = 1, <Vj' £ ) B > 0. Since 9 * (1, . . . , 1) then there exists / e {1, . . .,n} such 
that 

u(x + ei) = 6iu{x) for x = (xi, x 2 , . . . , - 1/2, . . . , x„), where 0/ * 1 (84) 

T 

/-th place 

We denote by S f the faces of Y which are orthogonal to the axis x t that is 

St = |xe5F:x / = ± i| 

We need an additional estimate 
Lemma 2.4. For any v 6 H l (F) the following inequality holds: 

\(v) sf -(v) F \ 2 <C\\Vv\\ 2 L2(F) (85) 



(here (v) s ± = \Sf\ 1 J v(x)dS x , where dS x is the volume (area) form on Sf, \Sf\ = J dS x ). 

sf sf 

Proof. Let v be an arbitrary function from C°°(F). One has 

a 

Cdv 

v(x) - v(x + ae/) = - — (x + e t t)dt (86) 
o 

where x e Sj, < a < (l-b)/2 (and therefore x+aei 6 F). We denote F = {x e F : -1/2 < x„ < -b/2}. 
Integrating equality (|86l ) by x\, x 2 , x/_i, x/+i, jc„ from -1/2 to 1/2 and by a from to(l-Z?)/2, 
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dividing by (1 -b)/2 and squaring we get 

|<v> 5r - <v) ? | 2 = 



1/2 1/2 (l-2>)/2/ a 

I -J I I d i (x+e " )d ' 



1/2 



-1/2 VO 



dadxi . . . dxi-idx, 



i+i 



< 



< C\\Vv\\l 2(F) (87) 



The fulfilment of inequality (f8~7l) for v e H l (F) follows from the standard embedding and trace 
theorems. 

It is well-known (see, e.g, Il23l Chapter 4]) that the operator II : H l (F) —> H l (Y) exists such that 
for an arbitrary v e H l (F) one has 



nv| F = v, ||nv|| ff i (y) < c||v|| ff i (F) 

We have the estimate which follows directly from ffT8l Lemma 3.1] 



|(v) ? -<v) f | 2 <ciivm 



\L 2 (Y) 



(88) 



(89) 



Then inequality ([85]) (with Sj) follows from (|8"71> -(l8"9l. The proof of ([83) with 5+ is similar. 
Lemma |24l is proved. □ 

Now using Lemma [241 and (f84l) we get 

|<vf >, | 2 = |i - ef 2 |<vf ) F - e^hf < 

< 2 |1 - 9,r 2 (|<vf £ > f - <vf ) s+ | 2 + Q 2 |<vf ) s _ - <vf%| 2 ) < C||V 
It follows from (f90b and the Poincare inequality that 



e,ei,2 

1 ll^2(F) 



(90) 



v 



e, £ |,2 
i Hl 2 (F) 



Ivf - <vf>rllLcp) + 



|<vf) F | 2 • |F| < C||Vv 



1 Hii(F) 



(91) 



Thus similarly to the Dirichlet eigenfunction v t ' e (see (1291) ) the function Vj' £ satisfies the Friedrichs 
inequality in F (although v e { £ 4- on dY\). As for the rest the proof of (1831 repeats word-by- word 
the proof of (1201) : we approximate the eigenfunction Vj' £ by the function vf £ (1481 (since vf £ = 
in the vicinity of dY then it satisfies (fTTT)-(fT2l) with an arbitrary 6) and then prove the asymptotic 
equality 



j0,e 



IVv 



D,£||2 



' i NL 2 (y E )' $ ^ (1> • • • > 1) 

The proof of (1921 (taking into account (1911 ) is completely similar to the proof of (l60t . 
Lemma [231 is proved. 

2.5. End of the proof of TheoremED It follows from CUB and O that 

oo 

(r^) = [jK( £ ),a { + ( g )] 

where the compact intervals [cT k {s), a k (e)] axe defined by 

[a- k (s),at(£)] = {J{s- 2 A d k t 

9eT" 



(92) 



□ 



(93) 



(94) 
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We denote B x = (1, 1, . . . , 1), 6> 2 = -0 X . It follows from Q3) and d94]) that 

e- 2 ^' e < aj(s) < e" 2 ^ (95) 

£ - 2 i^ £ < atis) < s- 2 A D k < £ (96) 

Obviously if k = 1 then the left and right-hand- sides of (1931 ) are equal to zero. It follows from 
(1621 ), (1821) that in the case k = 2 they both converge to /j. as s — » 0, while if & > 3 they converge to 
infinity in view of (ED), dgD, dSD- Thus 

a\{s) = 0, lima2(e)=//, lima^(e) = <x>, = 3, 4, 5 . . . (97) 

£^0 £— >() 

Similarly in view of (US), (Ug), (JSD), (ED, €3 one has 

limaj"(e) = o-, lim at (e) = oo, = 2,3,4 .. . (98) 

£^0 £— »0 

Then ©-© follow directly from (|93>, (I97l-(|98T>. Theorem Q is proved. 

3. Application to the theory of 2D photonic crystals 

In this section we apply the results obtained above to the theory of 2D photonic crystals. 

Photonic crystal is a dielectric medium with periodic structure whose main property is that the 
electromagnetic waves of a certain frequency cannot propagate in it without attenuation. From 
the mathematical point of view it means that the corresponding Maxwell operator has gaps in its 
spectrum. We refer to [|4|-(6) for more details. 

It is known that if the crystal is periodic in two directions and homogeneous with respect to the 
third one (so-called 2D photonic crystals) then the analysis of the Maxwell operator reduces to 
the analysis of scalar elliptic operators. In the case when dielectric medium occupies the entire 
space the rigorous justification of this reduction was carried out in IfTOll (in this work spectral 
gaps open up due to a high contrast electric permittivity). In the current work we consider the 
dielectric medium with a periodic family of perfectly conducting trap-like screens embedded into 
it. It means that we should supplement the Maxwell equations by suitable boundary conditions 
on these screens. In this case the analysis of the Maxwell operator reduces to the analysis of the 
Neumann or Dirichlet Laplacians in a 2-dimensional domain which is a cross-section of the crystal 
along periodicity plane. 

In this work we derive this reduction on a formal level of rigour. Then using Theorem II .H and 
Lemma |3~T1 below we open up gaps in the spectrum of the Maxwell operator. 
Let us introduce the following sets in R 3 : 

Sf = {(x l ,x 2 ,z) : x = (x u x 2 ) 6 S^, z e R} , £l £ = {(x l ,x 2 ,z) : x = (x u x 2 ) e z e R} 

where S £ and Q, £ belong to R 2 and are defined by © and © correspondingly. We suppose that 
fl £ is occupied by a dielectric medium while the union of the screens is occupied by a perfectly 
conducting material. It is supposed that the electric permittivity and the magnetic permeability of 
the material occupying Q e are equal to 1 . Then the propagation of electro-magnetic waves in Q, E is 
described by the Maxwell equations 

curlE = -— , curlH = — , divH = divE = 

ot ot 

supplemented by the following boundary conditions on [J S s : 

ieZ 2 

E r = 0, H v = 
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Here E and H are the electric and magnetic fields, E r and H v are the tangential and normal com- 
ponents of E and H correspondingly. 

Looking for monochromatic waves E(x, t) = e ,0>t E(x), H(x; t) = e mt H(x) we obtain 

AfU = coXJ, divH = divE = in Q. E , E r = 0, H v = on [J Sf (99) 

where U = (E, H) r , M E is the Maxwell operator: 

kas i -cuxl 
M=l \cuA 

(the subscript e emphasizes that this operator acts on functions defined in Q. E ). For more precise 
definition of the Maxwell operator we refer to [1241 (the theory developed in this paper covers, in 
particular, domains with a "screen-like" boundary). 

We are interested on the waves propagated along the plane z = 0, i.e. when (E,H) depends on 
x = (x\,x 2 ) only: 

E = (Ex(x u x 2 ), E 2 (x l: x 2 ), E 3 (x u x 2 )) , H = (Hi(x u x 2 ), H 2 (x { , x 2 ), H 3 (x u x 2 )) (100) 
Also we suppose that U is a Bloch wave that is 

36 = ...,#„)£ T : U(x + i) = ffU(x), i E Z" (101) 

(here ff = Y\(0 k ) ik ). We call the set of a> for which there is U = (E, H) + satisfying (|99l- (|101|) a 

(Bloch) spectrum of the operator M E . We denote it by cr(M E ). 
Using (1 1001) we can easily rewrite the equality M £ \J = oAJ as 

dH 3 dH 3 (dH 2 dHA 
-z'-— = dE u i—— = coE 2 , -i — - = a>E 3 (102) 

OX 2 OX\ \ OX\ ox 2 ) 

.dEi Mi AdE 2 dEA 

i- — = u>H u -i-T— = <x>H 2 , i — - = a)H 3 (103) 

OX 2 OX\ \ OX\ ox 2 ) 

Let us show that 

CO E Cr(M E ) <=> CO 2 6 (r($Q U Cr(J[ E ) (104) 

where ${i and J{ £ are, correspondingly, the Dirichlet and Neumann Laplacians in Q. E . Suppose that 
co £ cr(M E ). If co = then co 2 6 cr(Jl E ) otherwise we express Hi and H 2 (resp. E\ and E 2 ) from 
the first two equalities in (11031) {resp. (11021) ) and plug them into the third equality in (11021) (resp. 
(11031) ). As a result we get the following equalities on Q. E : 

-AE 3 = co 2 E 3 (resp. - AH 3 = co 2 H 3 ) (105) 

Let n = (ni,n 2 ,0) be the unit normal to US^. Since E T = then E\\n and therefore E ± (0,0, 1) 
and E ± (-n 2 ,ni,0). Taking this and the first two equalities in (11021 ) into account we obtain the 
following boundary conditions for E 3 and H 3 on |J S E : 

ieZ n 

OH dH dH 

E 3 =0, —l = —ln l + —^n 2 = ico- l (-E 2 n l +A e 2 E ~p.s 2 if 6 = (1, 1, . . . , 1)) = (106) 
on ox i ox 2 

Furthermore it follows from (I1021 i- (ll03t that 

U ± <=> E 3 t or H 3 t (107) 
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It follows from (1105l) - (1107l) and (11011) that co 2 belongs to the spectrum of either the Dirichlet or 
Neumann Laplacian in Q. £ . The converse implication in (11041 ) is proved similarly. 

Remark 3.1. Above we have dealt with the space 

J = {(E, H) : divE = divH = on £l £ , E T = H p = on U Sf and (fTOOl) - (fTOTT) holds} 

We introduce the following subspaces 

J E = {(E, H) eJ:E 1 =E 2 =H 3 =0}, J H = {(E, H) £/: H x = H 2 = E 3 = 0} 

Their elements are usually called E- and //-polarized waves. These subspaces are L 2 -orthogonal 
and it is easy to see that each U e J can be represented in unique way as 

U = U E + Uh, where U E e J E , U H e J H 

Moreover J E and J H are invariant subspaces of M £ . Thus cr(M £ ) is a union of cr(M £ \j E ) (E- 
subspectrum) and cr(M £ \j H ) (//-sub spectrum). We have just shown that the set of squares of points 
from /t-subspectrum is just the spectrum of the Dirichlet Laplacian, while the the set of squares of 
points from //-subspectrum is the spectrum of the Neumann Laplacian. 

The spectrum of Jl £ has been studied above (Theorem ll.il) . Now, we describe the spectrum of 
the Dirichlet Laplacian in Q £ . We define it via a sesquilinear form r] £ which is defined by formula 
© and the definitional domain dom^) = H^Qf) (here H^(Q. £ ) is a closure in H l (Q. £ ) of the set 
of compactly supported in Q £ functions). The Diriclet Laplacian in Q £ (we denote it Ji^) is the 
operator which is generated by this form, i.e. (0 holds (with J{ £ w t]^ instead of 3\ £ , if). It turns out 
that its spectrum goes to infinity as s — » 0, namely the following lemma holds true. 

Lemma 3.1. One has: 

mm{A : A e cr(^K)} -> oo (108) 

Proof. For i e Z" we denote B £ = s{B + i), F £ = s(F + i). Let u e dom^. Since u = on S £ one 
has the following Friedrichs inequalities: 

Nl^tBp < Cs 2 \\Vu\\ 2 L2(Bf) , WuWl^ < Ce 2 \\Vu\\ 2 L2(Fr) 

Here C is independent of u, i and e. Summing up these equalities by i e Z" and then integrating by 
parts we get: 

(3§u, u) Mm > Cie- 2 || M ||^ (ns) (109) 

It follows from (11091 ) (see, e.g., fl25l Chapter 6]) that the ray (-oo, C\s~ 2 ) belongs to the resolvent 
set of J\q that imply (1 1081 ). The lemma is proved. □ 

Thus using Theorem II .![ Lemma I3TI and (11041) we conclude that the spectrum of the Maxwell 
operator M £ has just two gaps in large finite intervals when s is small enough; as s — > these gaps 
converge to the intervals (- V<x, - ^fp) and ( Vcr, ^fp). 
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